Strong magnetic field enhancement of spin triplet pairing arising from 
coexisting 2kp spin and 2kp charge fluctuations 
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We study the effect of the magnetic field (Zeeman splitting) on the triplet pairing. We show 
generally that the enhancement of spin triplet pairing mediated by coexisting 2kF spin and 2kF 
charge fluctuations can be much larger than in the case of triplet pairing mediated by ferromagnetic 
spin fluctuations. We propose that this may be related to the recent experiment for (TMTSF)2C104, 
in which a possibility of singlet to triplet pairing transition has been suggested. 
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Spin triplet superconductivity is one of the most fas- 
cinating unconventional superconducting state. The in- 
vestigation of the mechanism of such a pairing state has 
been a theoretical challenge. Spin triplet pairing medi- 
ated by ferromagnetic spin fluctuations has been studied 
from the early days in the context of superfluid "^He, but 
another possibility has arisen for the past several years: 
triplet pairing mediated by coexisting 2k p spin and 2k p 
charge (or orbital) fluctuations proposed for Sr2Ru04[lj 
and for organic superconductors (TMTSF)2X (X=PF6, 
CIO4, etc.) BSSIiJ, 7] and 0-(BEDT-TTF)2l3. i 

For (TMTSF)2X in particular, the possibility of spin 
triplet pairing has been pointed out starting from the 



early days M, El, E, El, Q El, El, [H El, E 

and two of the present authors as well as several other 
groups have come up with the possibility of a close com- 
petition between singlet d-wave-like pairing and a triplet 
/-wave-like pairing due to coexisting 2kp spin and 2k p 
charge fluctuations. Very recently, an NMR study on 
(TMTSF)2C104 has pointed out a possibility of a tran- 
sition from a spin singlet pairing at low magnetic fields 
to triplet pairing or an FFLO state at high fields [2^1 . In 
fact, such a possibility of singlet to triplet pairing transi- 
tion under high magnetic field was pointed out theoreti- 
cally [Iiilll,!!!]. 

In the present Letter, we study the effect of the mag- 
netic field (Zeeman splitting) on the triplet pairing using 
random phase approximation (RPA). We show generally 
that the enhancement of spin triplet pairing mediated 
by coexisting 2k p spin and 2k p charge fiuctuations can 
be much larger than in the case of triplet pairing me- 
diated by ferromagnetic spin fiuctuations. Applying the 
idea to a microscopic model for (TMTSF)2X in which 
strong 2kp spin and 2fcj7 charge fluctuations take place, 
we actually show that the magnetic field enhancement 
of spin triplet /-wave pairing is strong compared to the 
enhancement of ferromagnetic-spin-fiuctuation mediated 
triplet pairing that occurs in a triangular lattice Hubbard 
model {24 . |25| . Due to this strong effect, we show that 
even when spin singlet pairing dominates in the absence 



of the magnetic field, a transition to triplet pairing may 
take place by applying the magnetic field. This is con- 
sistent with the possibility of the magnetic field induced 
singlet-triplet transition in (TMTSF)2C104[2l|. More- 
over, this strong magnetic field effect may be used as a 
general probe for identifying the pairing mechanism of 
triplet superconductors. 

The extended Hubbard Hamiltonian that takes 
into account the Zeeman effect is given as H = 

g^sB^.j^ ^sgn{a)(^^^Cicr , where creates an electron 
with spin a at the z-th site, tij represents the hopping 
parameters, U is the on-site interaction and Vij are the 
off-site interactions. q^ibB is the Zeeman energy with the 
spin quantization axis z || B. We ignore the orbital effect 
of the magnetic field. Within RPA, the effective pairing 
interactions mediated by spin and charge fluctuations are 
given as 



V' (q) 



{U + 2V{(i)Y 



Xch (q) , 



for spin singlet pairing, 

y*--(q) = T/(q)-2((7 + \/(q))y(q)x"' 



(1) 



(2) 



for spin triplet pairing with d 
= +f (-1) for a =T (i), 



z which means total 



(q) 



^(q) + ^xfp(q)-f/'xi-(q) 



{U + 2V{ci)r 



Xch (q) , 



(3) 



for spin triplet pairing with d || z {Sz — 0), where V (q) 
is the Fourier transform of the off-site interactions. The 
longitudinal spin susceptibility and the charge suscepti- 
bility are obtained by xfp — (x^^ + X^^ ^ X^^ ^ X^^) /2 
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and xch = (x^^ + X^^ + X^^ + X^O Z^- Here, 



-(4) 



and 



-xZ''{u + v^)xr 



(1 + Xr^q) (1 + Xr^q) -{U + Vq)- X^X^ 



-X5) 



where Vq stands for V (q) . The fongitudinal bare suscep- 
tibiUty is given as 

v^-rn^- (k + q))-/(e.(k)) 

^ e.(k + q)-e.(k) ' 

where (k) is the single electron dispersion that con- 
siders the Zeeman effect measured from the chemical po- 
tential /i, and /(^cr(k)) is the Fermi distribution func- 
tion. The transverse spin susceptibility is given as x — 
— Uxt~)^ where the transverse bare susceptibil- 
ity is obtained by 



/(C.(k + q))-/fe(k)) 

^,(k+q)-e^(k) ■ 



(7) 



To obtain the superconducting state, we solve the lin- 
earized BCS gap equation within weak-coupling theory. 



A'^(/)^(k) = -^^''(k-k') 



(kO)-/(g., (kQ) 

^,(k') + e.' (kO 



cj,^ (k'). (8) 



We consider singlet and triplet pairings with d || z 
(/i = s,taa) for opposite spin pairing (cr ^ a'), and 
triplet pairing with d ± z (/i = taa) for parallel spin 
pairing((T =t, |). (/>'^ (k) is the gap function and the 
critical temperature Tc is determined as the temperature 
where the eigenvalue A reaches unity. To give a refer- 
ence for the values of the magnetic field, we calculate the 
Pauh hmit by hbBp = UbkeT^V^- Although RPA 
may be considered as quantitatively insufficient for dis- 
cussing the absolute value of Tc, we expect this approach 
to be valid for studying the competition between differ- 
ent pairing symmetries. In fact, as we shall see, we find 
very good agreement between the RPA results and the 
already known results obtained by dynamical cluster ap- 
proximation (DC A) [25]. 

Before giving the calculation results, we show gener- 
ally using the above formula that the effect of the Zee- 
man splitting on the triplet pairing caused by the coex- 
istence of 2k F spin and 2kF charge fluctuations can be 
very special. First, let us consider a case where off-site 
repulsions are not present, so that only the spin fluctua- 
tions are relevant, and therefore possibility of triplet pair- 
ing superconductivity arises due to ferromagnetic spin 



fluctuations. In this case, the triplet pairing interac- 
tions reduce to V*""^ (q) = -U^x'^'^ and T^*'^* (q) = 

+^Xsp (q) ~ U'^xtp ("^t)' where the formula for x"" also 

reduces to 



X"" = 



(Tcr 

Xo 



i-f/'xrxr' 



(9) 



Here, we assume without losing generality that Xo" is 
enhanced while Xo'^ is suppressed by the magnetic field. 
(Whether a or a =| depends on the band structure 
and the band filling of the system as we shall see later.) In 
the first order of the magnetic field, xVp and xjp~ are not 
affected because exchanging t and | do not affect these 
quantities, x"^"^ is enhanced because the numerator Xq'^ 
in eq.® is enhanced, but again the term in the denomi- 
nator, Xo'^Xo'^j is not affected by the magnetic field in its 
first order. Thus, although y*'^'^ should dominate over 
its enhancement due to the Zeeman splitting oc- 
curs only through the direct enhancement of Xo i which 
may not be so large for realistic magnetic fields. 

When a possibility of triplet pairing arises due to the 
coexistence of 2kF spin and 2kF charge fiuctuations, 
where the latter are induced by off-site repulsions, the 
situation can change drastically. To make the discussion 
simple, let us consider a case with —{U+2V{Q,2kF)) — U, 
or equivalently U + V{Q2kF) — 0, for which, in the 
absence of magnetic field, xfp(Q2feF) = xi~(Q2feF) ^ 
Xc/.(Q2fc^) and thus F*'^'^(Q2fe^) = V'^^Q2kF) - 
— V'(Q2fep)- (We shall see later that our idea works for 
more general cases). Here, the singlet and the triplet 
pairing interactions have nearly the same absolute val- 
ues because the 2kF spin and the 2kF charge contri- 
butions work constructively (destructively) in the spin 
triplet (singlet) pairing interaction. [l|, 0, 0, 0, H, 0|- In 
this case, x'^'^ at q = Q2kF *^an be given by a reduced 
form 



x""(Q2fc,) 



xr 



(Q2feJ 



i + i^(Q2fcJxr(Q2/cp) 



(10) 



Here again, we assume without losing generality that 
Xo'^(Q2feF) (Xo'^(Q2feF)) is enhanced (suppressed) by 
the Zeeman splitting, xlp and xtjT again are not af- 
fected by the magnetic field in its first order, so V 
and V*°''^ are unaffected, while x*^*^ is again affected. 
The difference from the case with ferromagnetic spin 
fiuctuations lies in that since V{Q2kF) < 0; the de- 
nominator of eg. pH)) decreases as Xo'^ increases. The 
enhancement of x"^"^ due to this effect can be very 
large in the vicinity of 2fci? CDW ordering because 
1 + {U + 2y(Q2fe^))xo(Q2fep) = signals this order- 
ing, which is the same as 1 -I- T^(Q2feF)xo(Q2fcF) = 
when U + V^(Q2fcj?) = 0. Therefore, when the possibihty 
of triplet pairing arises in the vicinity of coexisting 2k f 
CDW and 2kF SDW phases, triplet pairing with parallel 
spins can be strongly favored by the Zeeman splitting. 
In actual cases, superconductivity is usually degraded 



3 



by the orbital effect under magnetic fields, but even in 
that case, the enhancement of triplet pairing due to the 
above effect should make the suppression to be moder- 
ate. For quasi-lD systems in particular, where the addi- 
tional node in the SC gap required in the triplet pairing 
does not intersect the Fermi surface (see Figlljb)), the 
coexistence of 2kF spin and 2kp charge fluctuations al- 
ready results in a subtle competition between singlet and 
triplet pairings 0, 0] , so that the strong enhancement of 
the triplet pairing interaction by the magnetic field may 
easily result in a singlet to triplet pairing transition. 
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FIG. 1: (color online) (a) The model on the triangular lat- 
tice (left), the /-wave gap function (center), and the density 
of states in the presence of the Zeeman splitting schemati- 
cally, the Fermi level is for a dilute band filling (right), (b) 
The model for (TMTSF)2X on a QID lattice (left), and d- 
wave (center) and /-wave (right) gap functions along with 
the Fermi surface The arrows represent the nesting vector of 
the Fermi surface. 

We now apply the above idea to actual systems. First, 
we consider a case where the possibility of triplet pairing 
occurs due to ferromagnetic spin fluctuations induced by 
the on-site repulsion. As a typical example, we consider 
a case on a triangular lattice with dilute band filling as 
shown in Fig. [TJa). In this case, possibility of spin triplet 
/-wave pairing due to ferromagnetic spin fluctuations has 
been pointed out previously 24 , 2^ . 

The band dispersion is given as ^o- (k) = 2t cos kx + 
2t cos ky + 2tcos{kx + ky) — fJ. + gi.iBBsgii{a). We take 
the transfer energy as the unit of the energy, i.e., t = 1.0 
The on-site interaction is U — 3.0, and the band filling is 
taken as n = 0.2. We take 128 x 128 fc-point meshes in 
the RPA calculation. When the magnetic field is absent, 
we obtain ksTc ~ 0.014. The Pauli limit corresponding 
to this Tc is ijlbBp ~ 0.017, which should be considered 
as a reference for the values of the magnetic field. When 
the Zeeman splitting is introduced, Xo^ becomes slightly 
larger than Xo because of the the increase of the density 
of states (DOS) at the Fermi level of the down-spin due 
to the Zeeman splitting (see FigHJ^a)). Thus, this corre- 
sponds to the case with cr =| in our general argument 
for the case with ferromagnetic fluctuations, so that we 



should focus on the enhancement of x^'^ ■ In order to 
measure how the magnetic field enhancement of is 
refiected to the enhancement of x^^i we introduce a pa- 
rameter aa{(\,B) = xr(q^^)/xrfa°) ' ^^^'^ 
^bB = 0.03 is plotted as a function of q. ~ 1 means 
that the effect of the denominator in cq. ([9]) is small, as ex- 
pected from our argument above. Effect of the magnetic 
field on the strength of the triplet pairing is shown in 
FigEJa), where we plot the eigenvalues of the linearized 
gap equation for each triplet pairings. As expected, due 
to the enhancement in ^ triplet /^^-wave dominates 
in the presence of the magnetic field. We note that this 
result for A closely resembles the field dependence of the 
pairing susceptibility calculated by DCA for the same 
system 25], suggesting the reliability of the present ap- 
proach. 
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FIG. 2: (color online) a^(q,B) at T = 0.02 for (a)the tri- 
angular lattice with n = 0.2 and (b)the QID model with 
V2 + Vpl = 0.85 = U/2. 



(a) 



1 

0.9 
0.8 
0.7 
0.6 



ST/ 



tt 



ST/' 



(b) 



ST/ 

0.01 0.02 0.03 

\iBB 



0.7 
0.6 
0.5 
0.4 
0.3 















SSd 






0.01 0.02 0.03 
\1bB 

ST/^^ in (b) 



ST/^''' in (a) 



0.01 0.02 0.03 
\1bB 

FIG. 3: (color online) The B dependence of A at T = 0.02 for 
(a)the triangular lattice for n = 0.2, (b)the QID model for 
V2 -I- Vpl = U/2. (c) The B dependence of the most dominant 
A for the cases shown in (a) and (b) normalized by A at B = 0. 

We now turn to the case where the possibility of 
triplet pairing arises due to the coexistence of 2kF spin 
and 2kF charge fluctuations. As a typical example, 
we consider the case of (TMTSF)2X. % i, d, i, 
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We adopt a 3/4-filled QID extended Hubbard model as 
shown in Fig[T]Jb).[^ The band dispersion is given by 
^cr (k) = 2^2; cos kx + 2ty cos ky — ^ + gfiBBsgn{a). We 
take t.j; = 1.0 as the unit of the energy, and ty — 0.2. 
As for the interaction parameters, we consider not only 
the on-site interaction U but also off-site interactions : 
nearest-neighbor(n.n.) Vi, 2nd n.n. V2, 3rd n.n. interac- 
tion V3 in the x-direction and n.n. interaction Vpi in the 
y-direction, where the Fourier transformed off-site inter- 
action is given as V (q) = 2Vi cos {k^) + 2V2 cos {2kx) + 
2V3 cos {3kx) + 2Vpi cos (ky). We set the on-site and the 
off-site interactions as U — 1.7, Vi — 0.9 and V3 = 0.1, 
and V2 and Vpi are varied. The band filling is taken as 
n = 1.5 in accord with (TMTSF)2X. We take 256 x 128 
k-point meshes in the actual calculation. We first con- 
sider the case with V2 = 0.45, Vpi = 0.4, for which 
U + V{Q2k,.) = (Q2fe^ ~ (7r/2,7r)) and thus singlet 
d-wave and triplet /-wave pairings are nearly degenerate 
in the absence of the magnetic field. As discussed in the 
previous studies [1, H, S I3| j the degeneracy is because the 
triplet and singlet pairing interactions are nearly equal at 
Q2fep , and the additional gap node in the /-wave pairing 
does not intersect the Fermi surface, so that the nodal 
structure on the Fermi surface is the same between d 
and /, as shown in Fig[ljb). Introduction of the Zeeman 
splitting enhances Xo^(Q2feF) suppresses XoHQ^kp) 
because the up-spin band becomes close to half filling and 
the electron-hole symmetry is more restored. Thus, this 
corresponds to a =t in our general argument for the case 
with coexisting 2kF spin and 2kF charge fluctuations, so 
we should look at the enhancement of due to the 
magnetic field. In Fig[2I^b), is plotted at hbB — 0.03, 
which largely exceeds unity at Cl2kp as expected from 
our previous argument that the denominator in eq. (|10p 
becomes close to 0. In Figl3l[b), we show the magnetic 
field dependence of the eigenvalue A. It can be seen that 
triplet /^^-wave is strongly enhanced due to the strong 
enhancement of X^^(Q2feF)- In Figl3l[c), we compare the 
B dependence of A normalized by its value at _B = for 
the two cases. We can see that the enhancement of the 
triplet pairing mediated by coexisting 2kp spin and 2kp 
charge fluctuations is much larger. 

In the above, we considered a case where XchiQ2kF) = 
Xsp(Q2feF) and thus the singlet d-wave and triplet /- 
wave pairings are nearly degenerate a,t B ~ 0, but even 
when singlet pairing dominates at B = a transition 
to triplet pairing can take place within realistic values 
of B due to this strong enhancement of the triplet pair- 
ing interaction. To see this in detail, we consider a case 
with V2 = 0.4, Vpi — 0.4, for which c?-wave dominates 
over /-wave for B = 0. Here, Tc = 0.012, which corre- 
sponds to the Pauli limit of jisBp ~ 0.015. We obtain in 
Fig. Ua) a pairing "phase diagram" in the ( V2 -I- Vpi )-B 
space obtained by comparing A at ksT = 0.012. The 
phase diagram for the superconducting state (spin singlet 
d-wave and triplet /^^-wave) and the normal state in the 
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FIG. 4: (color online) (a) Pairing phase diagram at fcsT — 
0.012. (b) Calculated T-B phase diagram in V2 = 0.4 and 
Vpi = 0.4. 



temperature-magnetic field space is shown in Fig. [H^b). 
Applying the magnetic field, the spin singlet d-wave (in- 
dicated by SC-SSd) gives way to the spin triplet /-wave 
pairing with S, = +1 (SC-ST/TT). Note that the Tc in 
the spin triplet /^^-wave channel increases with B be- 
cause we ignore the orbital effect. We expect that the 
orbital effect actually suppresses the Tc, but even in that 
case, the effect of the Zeeman splitting should strongly 
favor the occurrence of triplet pairing over singlet pair- 
ing. 

In conclusion, we have generally shown that the mag- 
netic field enhancement of the spin triplet pairing due to 
the coexistence of 2kp spin and 2kp charge fluctuations 
can be extremely large compared to that mediated by 
ferromagnetic spin fluctuations. Thus, even when spin 
singlet pairing dominates in the absence of the magnetic 
field, a transition to spin triplet pairing can take place by 
applying (not unrealistically large) magnetic field. This 
is consistent with the possibility of the magnetic field 
induced singlet-triplet transition in (TMTSF)2C104f2o'|. 
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